On (Not)-Constraining Heavy Asymmetric Bosonic Dark Matter 
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Recently, constraints on bosonic asymmetric dark matter iiave been imposed based on the ex- 
istence of old neutron stars excluding the dark matter masses in the range from ~ 2 keV up to 
several GeV. The constraints are based on the star destruction scenario where the dark matter 
particles captured by the star collapse forming a black hole that eventually consumes the host star. 
In addition, there were claims in the literature that similar constraints can be obtained for dark 
matter masses heavier than a few TeV. Here we argue that it is not possible to extend to these 
constraints. We show that in the case of heavy dark matter, instead of forming a single large black 
hole that consumes the star, the collapsing dark matter particles form a series of small black holes 
that evaporate fast without leading to the destruction of the star. Thus, no constraints arise for 
bosonic asymmetric dark matter particles with masses of a few TeV or higher. 
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The last few years stellar observations have been used 
in order to constrain specific dark matter candidates or 
predict interesting phenomena related to dark matter 
|l3|. Specifically, compact stars such as white dwarfs 
and neutron stars have been found to impose severe con- 
straints on some dark matter models |14l427| . In princi- 
ple, there are two types of effects that can take place in 
compact stars and can give rise to constraints on dark 
matter. The first typ e is related to the thermal evolution 
of compact stars 0[H,[l^[23|- In this case, annihilation 
of trapped weakly interacting massive particles (WIMPs) 
inside a compact star can produce significant amount of 
heat that can change the thermal evolution of the star at 
later times. As a result, stars old enough to be quite cold 
might maintain higher temperature due to the released 
heat. 

The second type of constraints is related to asymmet- 
ric dark matter [li. [2l[23l427j . In this case WIMPs carry 
a conserved quantum number and there is an asymme- 
try between the populations of WIMPs and anti- WIMPs, 
so that the annihilation is impossible in the present-day 
universe where only the WIMPs remain. Such kind of 
WIMPs can accumulate in a compact star quite effi- 
ciently as long as the WIMP-nucleon cross section ex- 
ceeds a certain critical value which in the case of a neu- 
tron star is quite small {a ^ 10^*^ cm^) [l^, several 
orders of magnitude smaller than current limits from di- 
rect searches. Under certain circumstances, the amount 
of WIMPs accumulated during the lifetime of the star 
is sufficient to result in a gravitational collapse of the 
the WIMPs into a black hole. In the case of fermionic 



dark matter the amount of WIMPs needed for collapse is 
~ Mpj/m^, A/pi being the Planck mass and m the mass 
of the WIMP. Unless one assumes very heavy WIMPs 
and/or extremely high background dark matter density 
at the location of a neutron star, this amount of WIMPs 
is difficult to accumulate during the lifetime of the star 
(although the situation may change in the presence of 
attractive self- interactions [25|'). 

On the contrary, if WIMPs are of bosonic nature, there 
is no Fermi pressure and the amount of WIMPs needed 
for gravitational collapse is significantly lower [2^ , 
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This amount of WIMPs can be easily accumulated even 
by nearby known old neutron stars with the standard 
assumption about dark matter density near the Earth. 
As it was pointed out in [23], once ~ 10^'' WIMPs have 
been accreted and thcrmalized within the neutron star, 
a Bose-Einstein condensate (BEC) forms and all newly 
accreted WIMPs fall into the ground state. This state 
is very compact, so the WIMPs in the condensate start 
to self-gravitate way before the condition of Eq. ([Ij is 
satisfied. Thus, once the amount of WIMPs in the con- 
densate meets the criterion for the gravitational collapse 
(Eq. ([IJ), they form a black hole. Note that only a frac- 
tion of all WIMPs (namely, those in the BEC) goes into a 
black hole. If the resulting black hole starts to grow and 
can destroy the host star in a reasonable time, constraints 
on dark matter arise as follows from the well known ex- 
istence of old (more than a few billion years old) neutron 
stars. 

The fate of the black hole inside the neutron star is de- 
termined by the two competing processes: the accretion 
of surrounding nuclear matter that increases the mass of 
the black hole and scales as (Mbh being the black 
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hole mass), and Hawking radiation that scales as -Mg^. 
Therefore, there exists a critical value Man such that 
heavier black holes grow until the whole star is consumed, 
while lighter black holes evaporate completely. It follows 
from Eq. (H]) that the mass of the black hole formed by 
this mechanism is inversely proportional to the WIMP 
mass, so that any black hole formed from WIMPs heav- 
ier than a few GeV is lighter than M^^it and evaporates 
completely, in which case no constraints arise. 

In Refs. [H m, 113] it was suggested that for WIMPs 
heavier than a few TeV the above mechanism is modified 
and the constraints are recovered. The argument makes 
use of the observation that for sufficiently heavy WIMPs, 
the self-gravitation of the WIMP sphere sets in before 
the formation of BEG. At this point the total amount of 
accumulated WIMPs exceeds by many orders the amount 
of Eq. ([Ij required for the black hole formation, as well 
as the critical value Merit- So if one assumes that the 
collapse happens to the whole WIMP sphere at once the 
resulting black hole is heavy enough to grow and destroy 
the star. Hence the constraints reappear. 

In this paper we demonstrate that this is not what 
happens. Although it is true that for heavy WIMPs the 
self-gravitation sets in without BEG formation, the col- 
lapse of the self-gravitating WIMP sphere is hampered 
by the released gravitational energy and happens on a 
timescale set by the cooling of the WIMPs. As the WIMP 
sphere gradually shrinks, its density increases and BEG 
forms on a timescale much shorter than the cooling time. 
From this stage on, the original scenario is reproduced: 
the BEG region grows, becomes self-gravitating and col- 
lapses into a black hole as soon as Eq. Q is satisfied. 
The resulting black hole is too small to grow. Instead, 
it evaporates on a time scale much shorter than it takes 
for the rest of the WIMP sphere to collapse. The net 
result is the formation of tiny black holes one after the 
other which evaporate shortly after their birth without 
producing a heavy stable black hole. Thus no constraints 
on heavy WIMPs arise. 

Let us consider this argument more quantitatively. 
The WIMPs become self-gravitating when their density 
becomes larger than that of the nuclear matter in the 
core of the neutron star. This happens when the total 
mass of the WIMPs exceeds 
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n.^2.2cmf^W^)-^^^ (3) 
VlO^Ky VlOTeV/ ' ^ ' 

is the thermal radius of the WIMP sphere inside the neu- 
tron star at temperature T and pc = 5 x 10^* GeV/cm^ 



is the star core density. The density of WIMPs required 
to form BEG is [Hi 

«^4.7.10»GeVc.„-(^)"^(^)"'. 

(4) 

which implies that the total WIMP mass necessary for 
the formation of BEG is 

MeEC^2.1xl„«GeV(j5!^^)(^)'. (5) 

From the dependence of Eqs. ^ and (O on the WIMP 
mass m one can see that for to > 10 TeV the self- 
gravitation sets in prior to the formation of BEG. Note 
that the total mass of the WIMPs in this case is large 
enough to overcome the constraint imposed by the un- 
certainty principle of Eq. ([iJ . 

Once the self self-gravitation sets in, the WIMP sphere 
starts to collapse. However, it cannot collapse directly 
into a black hole because the WIMPs have to lose their 
energy and momentum. Instead, the collapse of WIMPs 
is likely to resemble the formation of dark matter halos, 
with the essential difference that constant interactions 
with nucleons provide an extra energy loss mechanism 
for the WIMPs. One may expect, therefore, that the 
WIMPs develop a cuspy profile similar to the dark matter 
halos, which shrinks as the WIMPs lose their energy in 
interactions with nucleons. 

Even without the cusp, the shrinking WIMP sphere 
would form BEG long before it reaches the size compa- 
rable to its Schwarzschild radius. To make the argu- 
ment as robust as possible, consider the worst (unreal- 
istic) case where the WIMP sphere contracts maintain- 
ing a uniform density. In order to form a black hole of 
mass M the WIMPs have to reach the density pbh ^ 
3(327rG3M2)-i - 10^^ GeV/cm^(7\f/104°GeV)-2, while 
the density required for BEG formation is much lower, 
c.f. Eq. Note that the smaller is the black hole, the 
higher density is required, so the uniform contraction is 
indeed the worst case. Thus, BEG will unavoidably be 
formed before the WIMPs collapse into a black hole. 

Gonsider now the formation of BEG in some more de- 
tail. Once the density of Eq. ^ is reached, the formation 
of BEG happens on time scales of order [2^ 

iBEC -T^- lO-^'s, (6) 

i.e. practically instantaneously. Further shrinking of 
the WIMP sphere results in increasing the mass of the 
condensate rather than the density of non-condensed 
WIMPs. 

The size of the condensate can be estimated by not- 
ing that once the WIMPs become self-gravitating they 
dominate the gravitational potential. Taking their den- 
sity approximately constant and equal to that required 



3 



for the condensate formation gives 
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-3/8 
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This is much smaller than the size of the WIMP sphere. 

As the mass of the BEC grows, it eventually becomes 
self-gravitating itself. This happens when the density 
of the condensate becomes equal to the non-condensed 
WIMP density, i.e., when its mass exceeds 



Mbec, sg = yPcritriEC - 9.6 x lO^iGeV [j^) 

(8) 

Although self-gravitating, the condensate cannot yet col- 
lapse because it does not satisfy the constraint of Eq. ([IJ . 
The latter is satisfied when the total mass of the con- 
densed state reaches Mmin- Beyond that point stable 
configurations of the self-gravitating condensate do not 
exist and it collapses into a black hole [IHl • The resulting 
mass of the black hole is given by Eq. ([IJ . Note that this 
mass is much smaller than the total mass of the WIMP- 
sphere which is of order Msg , Eq. ([2]) . Note also that the 
black hole mass becomes smaller as the mass of the dark 
matter particle, m, increases. 

The black hole of mass Mj^in is too small to survive 
the Hawking radiation. It evaporated on the time scale 
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To conclude the argument, let us show that such burn- 
ing of dark matter by the formation and evaporation of 
small black holes is more efRcient than both the accretion 
of new dark matter onto the neutron star and the creation 
of new black holes by the above mechanism. To this end 
consider the relevant time scales. The accumulation of 
WIMPs by the neutron star proceeds at a constant rate 
F that depends on the local dark matter density pdm and 
equals 
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m 



F = 1.25 X lO^^GeV/s 
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103GeV/cm= 
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The time needed to accumulate the amount of dark mat- 
ter equal to Mmin is 



Consider now the time needed to form a new black hole 
in the course of collapse of the self-gravitating WIMP- 
sphere. According to the picture described above, the 
contracting WIMP-sphere reaches at some point the den- 
sity pcrit at which the formation of the BEC becomes pos- 
sible. Further contraction happens at a constant WIMP 
density p = Pcrit as the excess of the WIMPs goes into the 
BEC state. To form a black hole, the amount of WIMPs 
equal M„iin has to be moved from the WIMP-sphere to 
the BEC. The corresponding release of potential energy 
is equal to the black hole mass times the difference of 
the gravitational potential between the surface and the 
center of the sphere. 
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where M is the total WIMP mass and rg is the radius of 
the WIMP-sphere, 
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The released energy has to be dispersed via collisions of 
WIMPs with nucleons, which sets the time scale of the 
process. Taking into account the Pauli blocking factor, 
the WIMP-nucleon collision time is 



Tcol = 
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where pp is the Fermi momentum of the nucleons, niN — 
1 GeV is the nucleon mass, v is the WIMP velocity and 
e is the WIMP kinetic energy estimated as 



GMn 
2ro 



(15) 



from the virial theorem. Taking into account that the 
typical energy loss per collision is Se = 2mAre/m and 
assembling all the factors, the time required to lose the 
excessive energy and form a black hole is 
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One can see that this time is shorter for larger mass M. 
Substituting M = Msg into Eq. ([T6l) one gets 



tc 



1.5 X 10-^s 



7.6 X lO-^s 
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This is larger by an order of magnitude than the evapo- 
ration time even for m = 10 TeV and even considering a 
10^ GeV/cm^ local dark matter density. For nearby neu- 
tron stars where the dark matter constraints are more 
reliable, this will give 5 orders of magnitude higher for 
10 TeV WIMP mass. 
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This time is shorter by a factor of a few than the black 
hole evaporation time of Eq. (jQ]). Note, however, the 
strong dependence of both quantities on the WIMP mass 
TO. Already for masses m > 13 TeV the black hole evap- 
oration time becomes shorter. 
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Three comments are in order. The first comment is 
related to the black hole evaporation. The fate of a black 
hole of mass A/bh is dictated by the equation 
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dt 
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where the first term in the left hand side corresponds to 
the Bondi accretion, and the second term corresponds to 
the Hawking radiation. Here C — 4:TrXpcG^/cl, A be- 
ing a constant of order one and Cg the speed of sound. 
As we mentioned earlier, using Eq. ([l} one can find that 
for WIMP masses above roughly 10 GeV, Hawking ra- 
diation wins over accretion and the evaporation time is 
given by Eq. However, the newly formed black hole 
can also accrete from the dark matter population that, in 
the case examined here, can have densities significantly 
higher than the nuclear matter density pc, as it can be 
seen from Eq. Q . The accretion of non- interacting colli- 
sionless particles in the non-relativistic limit onto a black 
hole is given by (soj 
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where pdm is Pcrit of Eq- Oi a-nd Vrx, is the average 
WIMP velocity far away from the black hole. Using the 
virial theorem, we take t^oo = \/ GM/r, M being the to- 
tal mass of the WIMP sphere and r the radius of the 
WIMP sphere. It is understood that for the first black 
hole M = Msg, and r = tq. We have checked that for 
WIMP sphere masses ranging from Eq. ([T]) to Eq. ([2]) 
the Hawking radiation dominates overwhelmingly over 
the WIMP accretion despite the large WIMP density, 
and therefore the black hole evaporation time is given 
accurately by Eq. (jH]). 

The second comment is that the cooling time derived 
above refers to the formation of the first black hole 
of mass given by Eq. ([1]). Subsequent black holes re- 
quire progressively longer times. This is easily seen from 
Eq. because tcooi scales as 1/M^/'^ (recall that tq 
scales as M^^^ from Eq. ([13])). Thus, the more black 
holes have formed and evaporated, the smaller is the re- 
maining WIMP mass M and the longer the time needed 
to form the next one. 

The third comment has to do with the neutron star 
temperature. As one can see in Eq. ([TT]) , there is a strong 
dependence in T. So it is essential that 10^ K is a good 
estimate for the temperature of the neutron star at the 
time the star has accumulated Msg mass of WIMPs. Note 
that Msg depends on the temperature, cf. Eq. The 
neutron star does not always has the same temperature, 
but becomes colder with time due to different cooling 
mechanisms. The cooling curves (i.e. the temperature of 
the star as a function of the age) for typical neutron stars 
can be found, e.g. in [3l| and [11]. Knowing the cool- 
ing curves, Msg{T) and the WIMP accretion rate from 
Eq. ([TU]), one can easily estimate the time it takes for the 
star to accumulate a total WIMP mass of Msg. We found 



that in all cases of interest (i.e. for m ranging from 10 to 
1000 TeV and pdm from 0.3 to 10^ GeV/cm^), the accu- 
mulation time is at least ~ 1.5 x 10^ years. By that time 
the temperature of the star should be below 10^ K, which 
justifies the use of this temperature in our estimates. 

Finally, let us examine one more factor that can poten- 
tially affect our estimates. As the first black hole evapo- 
rates, it emits particles via Hawking radiation that reheat 
the core of the star. In principle such a deposit of heat 
can increase the temperature of nucleons and WIMPs 
and thus could potentially change our estimates. We 
show here that the increase in the temperature is small 
and the whole effect is negligible. 

The energy produced by the black hole evaporation, 
once transferred to nucleons, propagates according to the 
diffusion equation. 



(20) 



where u is the energy density due to the black hole evap- 
oration, D is the diffusion coefficient and Q is the energy 
injection rate. 

Since the black hole is tiny and the mean free 
path of the Hawking-emitted particles is short (^ 
10""'^^ cm assuming a typical nucleon cross section un ^ 
7rl0~^^ cm^), one may take Q to be a delta-function in 
space, Q = Q{t)5{x). We will also assume that the en- 
ergy is injected at a constant rate. As we will see below, 
the time scale associated with the heat conductivity of 
nuclear matter is extremely short, so a constant rate is 
a good approximation for most stages of the black hole 
evaporation. 

The diffusion coefficient D — k/cv can be expressed 
in terms of the specific heat capacity cy — 1.4 x 
10-9GeV^(T/105K) of the nuclear matter [32] and the 
thermal conductivity k. The latter can also be found in 
Ref. [13, lO^GeV^ 

Assuming that the energy release by the black hole 
starts at < = 0, the solution to the above diffusion equa- 
tion is 



u{r, t) 



Q 



-erfc 



(21) 



The time-dependent factor saturates to one as soon as 
r/y/iDt <C 1 and the solution of Eq. ([21]) becomes sta- 
tionary at corresponding distances. The characteristic 
time scale for reaching the stationary solution at dis- 
tances of order of the size of the WIMP cloud is 



tst^ 2.9 X lO^^s , 

AD V 



This time scale is much shorter than icooi for all cases 
of interest. The same time scale determines the onset of 
the equilibrium once the heat source is switched off (the 
black hole has evaporated). 

When the stationary solution is reached, the tempera- 
ture of the nuclear matter is given by 



inkr 



(23) 
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where T^o — 10^ K. The first term describes the effect of 
the evaporating black hole. 

At distances comparable to the size of the WIMP- 
sphere r ~ ro, the change in the temperature during the 
black hole evaporation is ~ 10 K (m/10 TeV)^, where we 
have estimated Q by dividing the black hole mass by its 
evaporation time. However, in cases where icooi is longer 
than the evaporation time, it is more appropriate to es- 
timate the value of Q by dividing the mass of the black 
hole by the icooi of Eq. ([T7]). In this case, we find an 
increase in the temperature at most 

/TO \-8/3 , ^ 

Thus, the effect of the black hole evaporation on the tem- 
perature in the center of the neutron star is completely 
negligible. 

I. CONCLUSIONS 

We demonstrated that it is not possible to obtain 
constraints on asymmetric bosonic dark matter with 
masses in the TeV range or higher from the collapse of 
WIMPs into black holes inside neutron stars and de- 



struction of the latter. Although the self-gravitation 
of the accumulated WIMPs starts before the BEC for- 
mation, we showed that the WIMP sphere unavoidably 
passes through the stage of the BEC formation as it con- 
tracts. As a result, black holes with masses much smaller 
than the total WIMP mass form and evaporate in times 
shorter than both the accretion time scale and the WIMP 
sphere cooling time scale. This means that the WIMP 
sphere does not collapse to a single large black hole which 
would grow and consume the host star, but instead col- 
lapses piece by piece in such a way that every black hole 
that forms evaporates before another gets formed. The 
overall effect is the (unobservable) heating of the neutron 
star, but not its destruction. Thus no constraints on the 
TeV WIMP mass range can be imposed. 

We would like to stress that the scenario we have stud- 
ied is the most conservative one. We assumed a uniform 
constant density for the WIMP sphere. In reality, one 
should expect a more cuspy profile towards the center of 
the WIMP sphere. As a result, psEC will be achieved 
earlier, and the cooling time for formation of subsequent 
black holes can increase substantially because the earlier 
the BEC forms, the smaller will be the energy loss of 
WIMPs colliding to nucleons due to smaller differences 
in the temperatures of the two species. 



J. Casanellas and I. Lopes, Astrophys. J. 705, 135 (2009) 
[arXiv:0909.1971 [astro-ph.CO]]. 

M. A. Perez-Garcia, J. Silk and J. R. Stone, Phys. 
Rev. Lett. 105, 141101 (2010) [arXiv:1007.T421l [astro- 
ph.CO]]. 

J. Casanellas and 1. Lopes, Mon. Not. Roy. Astron. Soc. 

410, 535 (2011) 'arXiv: 1008.0646 ' [astro-ph. CO]] . 

I. Lopes, J. Casanellas and D. Eugenio, Phys. Rev. D 83, 

063521 (2011) arXiv: 1102.2907 [astro-ph.CO]]. 

J. Casanellas and I. Lopes, Astrophys. J. 733, L51 (2011) 

[arXiA^lT0 4.5465 [astro-ph.SR]]. 

J. Casanellas, P. Pani, I. Lopes and V. Cardoso, Astro- 
phys. J. 745, 15 (2012) arXiv:1109.0249 [astro-ph.SR]]. 
M. A. Perez-Garcia and J. Silk, Phys. Lett. B 711, 6 
(2012) arXiv:1111.2275 [astro-ph.CO]]. 
M. A. Perez-Garcia, J. Silk and J. R. Stone, AIP Conf. 
Proc. 1441, 525 (2012) arXiv: 1108.5206 [astro-ph.CO]]. 
L. Brayeur and P. Tinyakov, Phys. Rev. Lett. 109, 
061301 (2012) [arXiv:111 1.3205 [astro-ph.CO]]. 
F. locco, M. Taoso, F. Leclercq and G. Meynet, Phys. 
Rev. Lett. 108, 061301 (2012) |arXiv:1201.5387 [astro- 
ph.SR]]. 

C. J. Horowitz, arXiv:1205.354T] [astro-ph.HE]. 

F. Capela, M. Pshirkov and P. Tinyakov, 

larXiv:1209.6"02T1 [astro-ph.CO]. 

I. Lopes and J. Silk, Astrophys. J. 757, 130 (2012) 
arXiv: 1209.3631 [astro-ph.SR]]. 

I. Goldman and S. Nussinov, Phys. Rev. D 40, 3221 
(1989). 

C. Kouvaris, Phys. Rev. D 77, 023006 (2008) 
,arXiv:0708.236 2 [astro-ph]]. 



[16 
[17 
[18 
[19 
[20 
[21 
[22; 
[23; 
[24 
[25; 
[26 
[27 
[28; 



G. Bertone an d M. Fair bairn, Phys. Rev. D 77, 043515 

(2008) arXi^0709.1485' [astro-ph]]. 

F. Sandin and P. Ciarcelluti, Astropart. Phys. 32, 278 

(2009) arXiv:0809.2942 [astro-ph]]. 

M. McCuUough and M. Fairbairn, Phys. Rev. D 81, 

083520 (2010) a rXiv:1001.2737 [hep-ph]]. 

C. Kouvaris and P. Tinyakov, Phys. Rev. D 82, 063531 

(2010) arXiv: 1004.0586 [astro-ph.GA]]. 

A. de LavaUaz and M. Fairbairn, Phys. Rev. D 81, 123521 

(2010) arXiv: 1004.0629 [astro-ph.GA]]. 

C. Kouvaris and P. Tinyakov, Phys. Rev. D 83, 083512 

(2011) arXiv: 1012.2039' [astro-ph.HE]]. 

P. Ciarcelluti and F. Sandin, Phys. Lett. B 695, 19 (2011) 
arXiv:1005.0857 [astro-ph.HE]]. 

C. Kouvaris and P. Tinyakov, Phys. Rev. Lett. 107, 

091301 (2011) arXiv:1104.0382 [astro-ph.CO]]. 

S. D. McDermott, H. -B. Yu and K. M. Zurek, Phys. Rev. 

D 85, 023519 (2012) [arXiv:1103.54 72 [hep-ph]]. 

C. Kouvaris, Phys. Rev. Lett. 108, 191301 (2012) 



Reno and 1. Sarcevic, 



[arXiv:1111.4364 l [astro-ph.CO]]. 
T. Guver, A. E. Erkoca, M. H. 
arXiv:1201.2400 [hep-ph]. 

Y. -z. Fan, R. -z. Yang and J. Chang, larXiv:1204.2564l 
[astro-ph.HE] . 

E. W. Mielke and F. E. Schunck, Nucl. Phys. B 564, 185 
(2000) gr-qc/0001061 . 
[29] H. T. C. Stoof, Phys. Rev. A 45, 8398 (1992). 
[30] Ya.B. Zeldovich and I.D. Novikov in Relativistic As- 
trophysics,vol. 1, The University of Chicago Press, 
Chicago, 1972 

[31] D. G. Yakovlev, O. Y. Gnedin, M. E. Gusakov, 



A. D. Kaminker, K. P. Leven fish and A. Y. Pot ekhin, Prakash, Astrophys. J. 425, 802 (1994). 

Nucl. Phys. A 752, 590 (2005) [astro-ph/0409751| . 

J. M. Lattimer, K. A. van Riper, M. Prakash and M. 



